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The infinite dimensional half-filled Hubbard model can be mapped exactly with no additional 
constraint onto a model of free fermions coupled in a Z 2 gauge-invariant manner to auxiliary Ising 
spins in a transverse field. In this slave-spin representation, the zero-temperature insulator-to- 
metal transition translates into spontaneous breaking of the local Z 2 gauge symmetry, which is not 
forbidden in infinite dimensions, thus endowing the Mott transition of an order parameter that is 
otherwise elusive in the original fermion representation. We demonstrate this interesting scenario 
by exactly solving the effective spin-fermion model by dynamical mean-field theory both at zero and 
at finite temperature. 

PACS numbers: 71.10.Fd, 71.30.+h, 05.30.Rt 


I. INTRODUCTION 


Let us consider the single-band model 


Sixty years after Sir Nevill Mott first envisaged it,® the 
interaction driven metal-to-insulator (Mott’s) transition 
remains an object of considerable efforts aiming at its 
full comprehension. The main difficulty resides in the 
inability of approximate methods based on independent 
particles to properly capture a transition that involves 
only part of electrons’ degrees of freedom: their charge. 
This is frustrating since the common attitude is trac¬ 
ing back complex many-body phenomena to independent 
particle models, an approach on which many successful 
techniques are based, e.g., the mean-field theory or the 
density functional theory within the local density approx¬ 
imation. 

A tool that has been often employed to disentangle 
charge from other electron quantum numbers, thus mak¬ 
ing Mott’s localization accessible already at the mean- 
field level, is to associate local charge configurations to 
novel fictitious degrees of freedom. This requires supple¬ 
mental local constraints in order to project the enlarged 
Hilbert space onto the physical one. 

A well-known implementation of this idea is the slave 
boson technique.^® The slave-boson mean-field theory is 
indeed able to describe the metal-insulator transition.® 
In this language, however, the Mott transition is asso¬ 
ciated with a local U( 1) gauge symmetry that is main¬ 
tained in the insulator but spontaneously broken in the 
metal, which is physically not possible. Restoring the 
symmetry to assess how robust are the mean-field results 
requires adding quantum fluctuations, which is not an 
easy task.®^ 7 

In a recently proposed alternative route the slave 
bosons are replaced by slave Ising variables P^ 11 The ad¬ 
vantage is that the size of the enlarged local Hilbert space 
is only a few times bigger than the physical one, not in¬ 
finitely as in the slave-boson approach. In addition, the 
continuous C/(l) local gauge symmetry is replaced by a 
discrete Z2 one. 
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where hi = c\ a c ia , and, for later convenience, we 
added a constant term. The Hamiltonian TL corre¬ 
sponds to the conventional Hubbard model when Ui = U, 
Vi, and to the Anderson (Wolff) impurity model when 
Ui = U for i = io while Ui = 0 for any i 7^ io, with the 
impurity sitting at ig. The partition function Z corre¬ 
sponding to the Hamiltonian in Eq. 0 can be shown 11 
to be equivalent to 

Z = Tr(e“^) =Tr(e- /m * JJPi), ( 2 ) 

i 

where the slave-spin Hamiltonian is 

= XX a i a l c L c j> -X^i’ ( 3 ) 

i,j (J i 


with the Pauli matrices cr“, a = x,y,z, describing auxil¬ 
iary Ising spins, and the projection operator 


with fi,; defined as [see Eq. 0 ] 

n l = i-2 (h i -i) 2 . (5) 


We observe that Hi has eigenvalue +1 when hi = 1 and 
— 1 when hi = 0 or 2; this operator thus corresponds to 
charge fluctuations away from single occupancy of site 
i. Therefore P, is actually a projector onto the physical 
Hilbert space defined by erf = Hi. As anticipated, the 
Hamiltonian "H* possesses a local Z2 gauge symmetry 

ofSi of, cL-^SiC^, (6) 
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where Si = ±1. 

The equivalence of partition functions in Eq. ([2]) was 
proved by equating the perturbative expansion in Ui of 
Z with that of Z * PA It was observed that the two terms 
that define the projector Pj in Eq. Q have a very distinct 
role in perturbation theory!- 11 Specifically, the first term 
1/2 accounts for all diagrams in the perturbation series 
where U, t is applied an even number of times, while the 
second term erf f\/2 takes care of all diagrams where Ui 
is instead applied an odd number of times. 

This simple observation has a remarkable outcome. If 
Z is known to be an even function of some Ui , then the 
term of I2j/2 in Eq. Q plays no role and can be dis¬ 
carded. This result can be easily proven even without 
resorting to the perturbation theory. If Z(Ui) = Z(—Ui), 
also Z^U) = Z^-Ut). Since a? H^U) ejf = H.(-Ui) 
and erf Pj af = 1 — Pj, because the partition function is 
invariant under a unitary transformation it follows that 



FIG. 1: (Color online) Sketch of the phase diagram of the 
paramagnetic half-filled Hubbard model on a Bethe latticed 
D = 2t is half the bandwidth. At U/D = 2.7, the first- 
order metal-insulator transition occurs at the temperature 
T fts 0.0037 D. 


Z * (Ui, Pj) 


Z.(-U i ,P i ) = Z.(U i ,l-P i ) (7) 


z * 



Pj + (1-Pj) ^ 


— Z t 



The fortunate event Z{Ui) = Z(—Ui) occurs, in par¬ 
ticular, at particle-hole symmetry in the Anderson impu¬ 
rity model. Moreover, through the well known dynamical 
mean-field theory (DMFT) mapping between an Ander¬ 
son impurity model and a Hubbard model on a lattice 
with infinite coordination number,® it follows that also 
the latter possesses this property at particle-hole sym¬ 
metry. In other words, the partitions functions of both 
models at the particle-hole symmetric point can be calcu¬ 
lated through the partition function of the corresponding 
spin-fermion Hamiltonian without any constraint. 

This mapping was exploited for the Anderson impu¬ 
rity model in Ref. 1131 where the equivalent spin-fermion 
Hamiltonian was exactly solve d by the numerical renor¬ 
malization group (NRGjP^^, unveiling the remarkable 
fact that the mean-field breaking of the local Z 2 gauge 
symmetry is not spurious: it actually occurs in the ex¬ 
act solution. In the same work it was also argued that 
a similar phenomenon might occur even in the Hubbard 
model on infinite-coordination lattices, where Elitzur’s 
theorem^ does not apply hence a local Z 2 gauge sym¬ 
metry can be spontaneously broken.® 

Our aim here is just to verify this conjecture and exam¬ 
ine all its consequences. As a byproduct, the numerically 
exact solution of the lattice model Eq. ([3]) provides for the 
first time the chance to assay how reliable is the mean- 
field approximation when applied to slave-spin models, 
at least in the case of lattices with infinite coordination 
number. 


II. THE MODEL AND ITS MEAN-FIELD 
SOLUTION 

The model we shall consider is the simple half-filled 
Hubbard model 


Hub 


-~7j Y, ( c lcj* + H -c-) 

v <i,j> ° 

1 ]> 


( 8 ) 


with nearest-neighbor hopping on a Bethe lattice in the 
limit of infinite coordination number, 2 —> 00 . We focus 
on the fully frustrated case where the antiferromagnetic 
order does not appear and there is a Mott transition from 
a paramagnetic metal to a paramagnetic insul ator. The 
phase diagram of the model is sketched in Fig. jjjbd The 
transition line between the metal and the Mott insulator 
is first order at any T > 0 with a critical endpoint at 
T c ~ 0.026 where D = 2t is half the bandwidth. 

As mentioned, the partition function of the Hubbard 
model ^Hub can be also calculated through 

^Hub=( 2 ) Tr(e-^- f ), (9) 


where N is the number of sites and 'H s -{ is the spin- 
fermion Hamiltonian, Eq. ([ 3 ]), with Ui = U,\/i , namely 


^-f = —7= Y. a i aX j kW + H - c ) 


<i,j> a 


V rr Z 


( 10 ) 


It follows that the spin-fermion model Eq. (10) must also 
have the same phase diagram shown in Fig. |jj 
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The Hamiltonian 'H s .{ is invariant under the N mu¬ 
tually commuting unitary transformations of f2j, Vi. 
Therefore we can divide the full Hilbert space into 2 N 
orthogonal subspaces identified by the sets of bits {n} = 
(ni, n 2 , ■ ■ ■, rii ,..., tin), with rii = 0,1, which contain 
states such that 


In turn, | fermions) is the ground state of 

^fermions = Yh { C \a C jcr + H.C.). (16) 

* <i,j> o’ 


Here 


of | tp{ n }) = (~l) ni I V’m)- (11) 

It follows that 

Z s . f = Tr(e-^- f ) = £ Tr w (e-^- f ) 

M 

= Y, ^s-f{n}, (12) 

M 

where Tr/„i stands for the trace within the invariant 
subspace {n}. All invariant subspaces are actually de¬ 
generate at particle-hole symmetry and in lattices with 
infinite coordination. Therefore Z s _ f{ ra } = Z*, V {n}, 
so that Z s _ f = 2 n Z*; a degeneracy that cancels out the 
prefactor in Eq. ([9]). 

In view of the above result, it is worth briefly discussing 
the meaning of gauge symmetry breaking. We observe 
that, within each subspace, the matrix elements of any 
gauge variant operator vanish. For instance, 

<^{„ } I Of I V-In}) = (-l) 2ni (V>M I Of of of n< I V> W > 

= I Of | 

hence (4>'r n y \ of \ il>{ n y) = 0. Therefore spontaneous 
gauge symmetry breaking can be revealed within each in¬ 
variant subspace only by the long-time behavior of gauge- 
invariant correlation functions, e.g. lim t ^. 00 (crf (t) erf (0)). 
Alternatively, one can add to the Hamiltonian a test field 
that explicitly breaks the symmetry, and study the re¬ 
sponse in the limit of vanishing field. 


A. Zero temperature 


t*=t (Ising | afaf |Ising), 


(17) 


does not depend on the specific nearest-neighbor bond as 
we have assumed translational symmetry. 

It follows that | fermions) is just the Fermi sea of a 
tight-binding Hamiltonian with constant hopping. This 
implies that J in Eq. (151 has the value J = 8f/37r inde¬ 
pendent of |Ising). 


The Ising model in Eq. (14) with constant J can be 


readily solved on the Bethe lattice, where the mean-field 
results becomes exact for z —> oo. In particular one finds 
that, for any site i 


I (Ising | < | Ising) | = 8 (U c -U) Ji¬ 


ll 2 

U: 


2 ’ 


(18) 


where U c = 8J. Therefore, for U < U c , the mean-field 
state has a finite order parameter that vanishes at U c 
and above it. In the original model, Eq. (101, the finite 


value of (erf) actually implies that the local Z 2 gauge 
symmetry Eq. © is broken. In fact, we observe that, 
given the variational wavefunction | T) in Eq. (13), the 
state 


i*>=n fa «<)">>> (19) 

i 

with nrii = 0,1, is also solution of the mean-field equa¬ 
tions with the same energy and is orthogonal to | \k). In 
fact, for i ^ j the average value 

0F| Hifij l-F) = (Tj Q f |4')(’F| % |4r) + 0(z" 2|< - J ' 1 ) 
~ O (2" 2|i_j| ), ( 20 ) 


We shall start our analysis with a mean-field decom¬ 
position. At zero-temperature, we assume a factorized 
variational state 


| \k) =| Ising) (g) | fermions). (13) 

The Ising wave-function thus corresponds to the ground 
state of the Ising model in a transverse field 


n lsing = -- Y ■/,, of of -^Y ( 14 ) 

<i,j> i 

where, assuming translational symmetry, 

Jij = (fermions | J i(j c^ a + H.c. | fermions) 


= Ji 


V < i,j > . 


(15) 


because the average of any vanishes on the half-filled 
Fermi sea. Since for any given i there are z r sites j at 
distance \i — j\ = r, one readily concludes that | \k) and 
| \E ,/ ) are indeed orthogonal for z —> oo. Therefore the 
mean-field lowest energy state is 2 N times degenerate, 
once more signaling the gauge symmetry breaking. Such 
degeneracy exactly compensates the entropic pre-factor 
in Eq. & thus providing a physically meaningful vari¬ 
ational estimate of the Hubbard model ground state en¬ 
ergy. We note that, on the insulating side where f* = 0, 
see Eq. ( p~7|) , the fermionic mean-field Hamiltonian in 
Eq. (16) vanishes, hence any wavefunction | fermions) 
solves the variational problem, leading to a degeneracy 
4^ that is only partly compensated by the pre-factor in 
Eq. 0. The net result is a residual entropy S = N In 2, 
which physically corresponds to the spin entropy of the 
hypothetical paramagnetic Mott insulator. 
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At first sight, the spontaneous symmetry breaking 
might look wrong by the Elitzur’s theorem. 16 In re¬ 
ality, as we mentioned, this theorem does not apply 
in infinite-coordination lattices, where such a symme¬ 
try is not impeded from breaking spontaneously. ^ The 
phase with broken Zi gauge symmetry actually corre¬ 
sponds to a metal phase in the Hubbard model ([8]), while 
the phase at U > U c to the Mott insulator. Thence 
U c = 32/37r D ~ 3.39 D is the critical value of the Mott 
transition, which slightly overestimates the actual value 
U c ~ 2.97 D, see Fig. |T| 

We mention that the zero-temperature mean-field so¬ 
lution is equivalent to the Gutzwiller approximation ap¬ 
plied to the Hamiltonian Eq. ^ EUSl an approximation 
that in fact provides exac t variational results in infinite- 
coordination lattices JHEH The inclusion of quantum fluc¬ 
tuations on top of the mean-field solution is equivalent 
to the time-dependent Gutzwiller approximation. 23 Fol¬ 
lowing Ref. till if we apply the spin-wave approximation 
to the Ising model by writing 


~ sin 0 (2-x 2 -p 2 ) 

— V2 cos 9xi, 

(21) 

a i - cos e (2 -x 2 -p 2 ) 

+ ypl sin Oxi, 

(22) 

erf ~ V2pi, 


(23) 


where Xi and pi are conjugate variables and 


cos 0 


U/U c if U< U D , 
1 if U > U c , 


(24) 


for the original Hubbard model Eq. ([8| the density ma¬ 
trix 


PG = V Pfermions V^ = V\ Pfermions V\ , (28) 

i 

where Pfermions is the Boltzmann-Gibbs distribution of 
free fermions with nearest-neighbor hopping at half¬ 
filling, and the linear operator 



with and real variational parameters such that 
<F 2 + = 1. We observe that 1 + ft, and 1 — 

project onto states where site i is singly occupied and 
empty/doubly-occupied, respectively. If we formally 
transform 


V\ —> Vi — 


V2\ 


1 + (_1)7U 0 . |fl) 


-^(i-(-ir^) m, (30) 


where erf |l) (-U-)) = ± lit (1))), we realize that Eq. (30) 


applied on an electronic wavefunction generates a spin- 
fermion state belonging to the symmetry invariant sub¬ 
space with {n} = (ni, n, 2 , • ■ ■, n,,.. ., njv), see Eq. (11). 


Replacing Vi with V % in Eq. (28) yields a variational 


density-matrix for the spin-fermion Hamiltonian Eq. (10) 
restricted to the invariant subspace {n}: 


the Hamiltonian (|14|) becomes quadratic 


Rising — Vq + ^ (. 


b 2 
2 z 


Xi x i ’ 

<M> 


(25) 


where Eq is the mean-field energy, a = U c /2 and b = 
U 2 /2U C for U < U c , while a = U/2 and b = U c /2 for 
U > Uc-^ 1 The spin-wave spectrum is limited within the 
energy window < w < um, where u>m = \J a(a + b) 
and 

w m = y/a{a-b). (26) 


In other words, the Ising excitations are gapped every¬ 
where except at the Mott transition where a = 6, close 
to which 


Wm 


Uc 

2 



1/2 


(27) 


B. Finite temperature 

The above mean-field approach can be extended at fi¬ 
nite temperature without embarking into delicate issues 
related to gauge redundancy by simply exploiting the 
aforementioned equivalence with the Gutzwiller varia¬ 
tional approachPB Within the latter scheme, one assumes 


PG P{ n} 


IP. 




(31) 


Since at half-filling Tr ^pfermions = 0, it follows 

that on the Bethe-lattice with z —> oo, 

Tr(p { „}of) = = m 

and 

Tr(p {D} a*cla*c jrT ) = 4$ 2 $ 2 

Tl' ( ^fermions cj a Cj a j . 


Both results are independent of {n}. This entails a renor¬ 
malized hopping t —» 4 <F 2 t = (l — m 2 ) t. 

Following the finite-temperature calculations of 
Refs. [2H and [5S] obtained within the Gutzwiller approx¬ 
imation, one has ultimately to minimize the free-energy 
density 


F(m) = — 


4T 


J de VD 2 - e 2 In (l + e~ p t 1 "” 1 ') e ) 

(32) 


7 tD Z J_ 

(m- 1) +T In2 
4 v ’ 

l-m 

V 2 2 


1 + m 1 + m 
- In- 
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where we recognise that the last term is just the entropy 
of a spin in a magnetic field B = Ttanh -1 m directed 
along z. The phase diagram of Eq. ( [32] ) was calculated 
by the saddle-point method in Ref. 1241 and agrees quali¬ 
tatively with the exact one shown in Fig. |Tj 


Since the variational free-energy functional, Eq. (32), is 
the same for any of the 2 N subspaces {n}, such a degen¬ 
eracy cancels the entropic pre-factor in Eq. |9]), bringing 
once again a free-energy that is well defined in the ther¬ 
modynamic limit. The above simple variational calcula¬ 
tion is however unable to assess whether the symmetry is 


preserved or spontaneously broken. In Sec. Ill B we shall 
instead perform a numerically exact calculation at finite 
temperature to answer such question. 


III. DMFT SOLUTION 

In order to assess the validity of the mean-field approx¬ 
imation, we now move to the exact solution of the model 
Eq. (101 by means of the DMFT using the NR.gP^ES. 


solve the impurity problem defined by 
’Himp = ^ e k cjctr c kr7 

k<T 


u 


+a* J2 M c L d * + H.c .) - - a z , (33) 


k a 

which describes a localized d-level that hybridizes with 
a bath of conduction electrons. The sign of the hy¬ 
bridization depends on an Ising degree of freedom u x 
whose dynamics is controlled by a field (7/4 in the 21- 
direction. Within DMFT, the bath parameters ek and 14 
are self-consistently determined. 12 Specifically, through 
the Hamiltonian Eq. (331 we calculate the single-particle 


Green’s function of the physical electron a x d a in Mat- 
subara frequencies: 

G(zcon) = - fdre^ (T T (a x (r) d„{r) **(0)4(0))), 

(34) 

and we require that the self-consistency is reached when 


V . 14 =t 2 G(ico n ). 

“ iu n - e k 

k 


(35) 


We shall focus hereafter on three dynamical quanti¬ 
ties: the Green’s function of the physical electron defined 
above, Eq. (341, the one of the auxiliary fermion d a , 


G d (iuj n ) = -J^dre^ (T t [d a {r ) 4(0))), (36) 

and finally the Ising spin Green’s function 

G a (in x ) = — j^dr e iaxT (!T t (* x (t) **( 0))), (37) 


where oj n and Q\ are fermionic and bosonic Matsub- 
ara frequencies, respectively. The corresponding spectral 
functions are defined by the analytic continuation to the 
real axis, for instance 


Mu) = -z 


Qm G(iuj n 


z0+), 


and similarly for G d (iio n ) and G ct (^a)j to which we shall 
associate A d (u>) and H CT (o;), respectively. 

The NRG calculation have been performed using dis¬ 
cretization parameter A = 2, with twist-averaging over 
N z = 8 interleaved discretization grids. The calcula¬ 
tions can be performed taking into account the spin 
SU{2) symmetry and the axial charge (isospin) SU(2) 
symmetry. The finite-temperature expectation values 
and dynamical quantities h ave b een computed by the 
full-density-matrix algorithm^ 27 ^and the self-energy has 
been determined using the approach with an auxil¬ 
iary Green’s function defined for the physical electron 

operator.^ 

In order to reveal the spontaneous symmetry breaking, 
we shall add to the Hamiltonian Eq. (33) a test field 


SH = —h x <j x , 


(38) 


and calculate the average of a x as h x —» 0. Moreover, we 
shall separately discuss the solution at zero and at finite 
temperature. 


A. T = 0 DMFT results 


The impurity model Eq. (331 without self-consistency 


and with a constant density of states of the bath was 
studied at zero temperature in Ref. [T31 It describes a 
dissipative two-level systempSl where the eigenstates of 
a x label the two levels and a z induces tunneling among 
them. The dissipative bath is actually sub-Ohmic so that 
the tunneling a z is irrelevant and the level localizes, i.e. 
(a x ) 7^ 0.^ This case was argued to correspond to the 
metallic phase of the Hubbard model Eq. ([8]). The issue 
is whether this is indeed true, namely if the localization 
revealed by ( a x ) 7^ 0 survives the self-consistency condi¬ 
tion Eq. (|35| . 

In Fig.[2]we show the spectral functions of the physical 
and auxiliary fermions, A(uj) and A d (u)), black and red 
curves on the left panel, respectively, as well as of the 
Ising spin A CT (o;), right panel, for different values of U. 
Several features are worth being highlighted. From the 
physical fermion spectral function A(uj) we can locate the 
Mott transition between U/D = 2.8 and 3, in full agree¬ 
ment with the direct studies of the Hubbard modelP^The 
manifestations of such a transition in the behavior of the 
auxiliary fermion and the slave spin are quite interest¬ 
ing. The auxiliary fermion has a spectral function A d [u j) 
resembling that of a resonant model that shrinks contin¬ 
uously approaching the Mott transition, above which its 
weight is concentrated in a <5-peak that is not visible in 
the figure. 
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FIG. 2: (Color online) Left panels: zero-temperature spec¬ 
tral functions of the physical and auxiliary fermions, A{ui) 
and Ad{tjo) shown as black and red curves, respectively, for 
a range of U. Right panels: corresponding Ising spin spec¬ 
tral functions A a (uj), which, being the imaginary part of a 
causal response function, is odd in frequency. For U/D = 3 
the system is in the insulating phase. 



FIG. 3: (Color online) Ising magnetization (a x ) as a function 
of an external test field h x in the metallic phase at U/D = 2.6 
and at zero temperature. We note the finite value as h x —> 0. 


The Ising spin shows even more remarkable dynami¬ 
cal features. In the metallic phase the spectral function 
A a (uj) displays two peaks at the position of the Hub¬ 
bard bands but also a linear component A a (uj) ~ ui close 
to zero frequency. In other words, unlike in the mean- 
field solution, the Ising spectral function extends down 
to zero frequency. The linear low-frequency behavior oc¬ 
curs in an interval that shrinks as the Mott transition is 
approached, and meanwhile the slope increases, just as 



FIG. 4: (Color online) Ising magnetization (a x ) as function 
of an external test field h x in the metallic phase at U/D = 2 
(top panel) and insulating one at U/D = 4 (bottom panel), 
for different temperatures. 


we would expect to observe in the imaginary part of the 
physical-electron charge response function. Above the 
transition, this low energy component disappears, and a 
gap opens in the spectrum. This is consistent with the 
picture of the Mott insulating state where the charge de¬ 
grees of freedom freeze on a high energy scale. 

In Fig. [3] we plot the Ising magnetization (a x ) as a 
function of an external test held h x , see Eq. (38), in the 
metallic phase at U/D = 2.6. We observe that ( a x ) 
approaches a finite value as h x —I 0, thus revealing the 
spontaneous breaking of the Z 2 gauge symmetry. In the 
Mott insulating phase (cr x ) instead vanishes for h x —> 0. 
It thus follows that (a x ) is a legitimate order parameter 
of the zero-temperature Mott transition in the slave-spin 
language, finite in the metal and zero in the insulator, in 
qualitative agreement with the mean-held calculation. 


B. T > 0 DMFT results 


The zero-temperature results demonstrate that the 
speculations based on the two-level system Hamiltonian 
Eq. (33) are reliable. We observe that while at T = 0 a 


two-level system coupled to a sub-ohmic bath localizes, 
at any T / 0 it does not, i.e. (cr x ) = 0, VT > 0P*1 The 
recovery of the gauge symmetry at hnite temperature is 
evident in Fig. [4] where we show the slave-spin magneti¬ 
zation ( a x ) vs. h x in the metal (upper panel) for differ¬ 
ent temperatures. In the insulator (bottom panel), the 
magnetization for h x —> 0 is practically zero at any tem¬ 
perature. In the metal, however, we observe the gradual 
onset of a hnite order parameter as T —> 0 before h x 0. 
This is further confirmed by the behavior of the zero-held 
slave-spin susceptibility, x = d(cr x )/dh x . In the metal 
phase, for hxed and sufficiently small h x the susceptibil¬ 
ity is increasing as a power-law T -1 in the range T > h x 
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before saturating for T <C h x . The saturated value of 
x(T = 0) diverges as h~ x . In the insulator, the suscep¬ 
tibility at fixed h x is weakly increasing with decreasing 
T and it saturates on a temperature scale of order band¬ 
width. The saturated value of \(T = 0) does not depend 
on h x in the insulator. 



FIG. 5: (Color online) Retarded Green’s function Cj) lir] (w) of 
the Ising operator a x at U/D = 2.7, i.e. in the metal phase 
not far from the transition, for different temperatures. We 
note the gradual emergence of a 5-peak at oj = 0 on decreasing 
temperature toward T = 0. 



cn/D 


The establishment of spontaneous symmetry breaking 
is also clear in Fig. [5j where we plot (at zero test field, 
h x = 0) the correlation function of the Ising spin, i.e. 

^2 e ~ PEn \{ n | o x | m)\ 2 S(uj- E m + E n ) 

n 

/ oo 

dte iut {a x (t)a x ( 0)), (39) 

-OO 

which is better suited to capture the low-frequency sin¬ 
gularity without the spectral weight cancellation due 
to broadening that affects more severely the (odd in 
frequency) Green’s function A a (oj). We observe that 
^spin( w ) gradually develops a 5(w)-peak as T decreases, 
which means that 

lim lim {a x (t)a x ( 0)) —> (a x ) 2 > 0, 

t —>oo T—>0 

hence the spontaneous symmetry breaking.^ 

In Fig. |6](a) we show the temperature evolution of the 
spectral functions A(lj), A d (uj) and A a {uj) at U = 2.7 D , 
when the system is in the metal phase at zero tempera¬ 
ture and, upon heating, crosses the first-order transition 
and turns into a Mott insulator, see Fig. [I] The close-ups 
to the low-frequency behavior are shown as insets. We 
observe that, across the phase transition, the density of 
states of the physical electron jumps down from a finite 
to an extremely small value, as shown in detail in Fig. [6] 
(right inset to top panel), which is what one also obtains 
directly from the Hubbard model Eq. ([8]). 


FIG. 6: (Color online) Spectral functions of the physical elec¬ 
tron (top panel), of the auxiliary fermion (middle panel), and 
of the Ising spin (bottom panel) at U = 2.7 D upon increasing 
the temperature across the first order metal-insulator tran¬ 
sition. The insets show the close-ups to the low-frequency 
behavior. The right-hand inset in the top panel shows the 
density of states at the chemical potential of the physical elec¬ 
trons at U = 2.7 D upon increasing the temperature. 


The novelty here is represented by the behavior of the 
auxiliary fermion and slave-spin spectral functions. We 
observe that A d (co) does not change qualitatively across 
the transition; it always looks like the spectral function 
of a resonant level model, whose width and height simply 
jump at the transition. 

More remarkable is the behavior of the slave-spin spec¬ 
tral function A a (ui). Until the system stays in the metal 
phase, T < 0.0037 D, H^w) ~ uat small w, with a slope 
that is practically temperature independent, see Fig. [6] 
bottom panel. Unlike at the zero-temperature Mott tran¬ 
sition, the interval in which A a (ui) ~ to increases with 
temperature until, above the first order transition, sud¬ 
denly disappears; all low energy spectral weight is trans¬ 
ferred at high-energy, see Fig. [6j and a gap opens. 

For T > T c ~ 0.026 D 1 i.e. above the critical end¬ 
point, see Fig. [l] there is no metal-insulator transition 
anymore upon increasing U ; just a smooth cross-over be¬ 
tween a bad metal (large resistivity, increasing with T) 
and a Mott insulator (large but finite resistivity, decreas¬ 
ing with T). There is no unique definition of the crossover 
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C. Gauge-symmetry breaking’s fate at T > 0 
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FIG. 7: (Color online) (a) Physical-electron spectral function 
intensity at u> = 0 and the slave-spin susceptibility d(a x )/dh x 
at constant temperature T > T c as functions of U. The circles 
indicate the inflection points on both curves that may serve 
as possible definitions of the cross-over between the bad metal 
and the Mott insulator, (b) Auxiliary-fermion and (c) slave- 
spin spectral functions across the metal-insulator cross-over. 


location, but it is empirically found that different choices, 
using either the spectral function, the resistivity, or the 
double occupancy, bring to almost coincident values (e.g., 
see Fig. 4 in Ref. US). We have found that in the slave- 
spin language the cross-over also reveals itself by an in¬ 
flection point in the slave-spin susceptibility d(a x ) /dh X: 
see Fig. The evolution of spectral functions across 
the cross-over is smooth although we find a notable dif¬ 
ference between auxiliary fermions and slave spins. The 
auxiliary-fermion spectrum changes only qualitatively: 
with increasing U the width of the peak narrows down 
and its height thus raises, see Fig.[7[b). Since the spin op¬ 
erators of auxiliary and physical fermions map onto each 
other, the above evolution simply represents the gradual 
formation of localised moments. On the contrary, the 
slave-spin spectrum, which describes charge degrees of 
freedom, shows a more pronounced variation: the linear 
slope of the low-frequency component drops markedly, 
see Fig. [7)(c), signalling a rapid cross-over from a gapless 
to a pseudo-gap regime. 


The previous results confirm the expectation that, as 
soon as T > 0, the Ising order parameter (of) = 0, in the 
sense that its expectation value in the presence of a test 
field h x vanishes when h x —> 0. To state this in a different 
manner, we have found that the 2 N degeneracy of each 
many-body excited eigenstate is not split linearly in the 
field strength h Xl unlike what happens for the ground 
state. 


We observe that the model Hamiltonian Eq. (101, 


which describes fermion and spin fields coupled in a 
gauge-invariant manner, lacks an explicit gauge field. 
The latter could be made explicit should we decou¬ 
ple the spin-fermion interaction through a Hubbard- 
Stratonovich transformation involving an auxiliary link- 
variable that would just play the role of the gauge-field. 
A question we may then ask ourselves is whether, even 
though the ’’matter” field, i.e. the Ising held erf, has van¬ 
ishing expectation value, still the gauge held has a hnite 
average value, as one would expect in a Z 2 gauge-Higgs 
theory!' 11 1:1 In order to address this question we add to 
the Hamiltonian Eq. © another source held 


sn s - { = -^= © ( 


C la C ja + H - C -)- 


(40) 


<i,j> 0 


At hnite h the gauge variant operator c\ a c^ a acquires a 
hnite expectation value, and, accordingly, also (of erf) be¬ 
comes hnite. Within the DMFT mapping onto an impu¬ 
rity model, the source term Eq. (401 amounts to adding 
to Eq. (f33l) the additional term 


which corresponds to hybridizing the impurity with 
another conduction channel /, and impose, besides 
Eq. ([351) , also 


V 2 

y* v f k 
“ - e/ k 


h 2 G d (iu n ). 


(42) 


The impurity model 'Him P +6TLi mp supplemented by the 
self-consistency conditions Eqs. (35) and (42) can still be 
solved by NRG, even though it now involves two chan¬ 
nels of conduction electrons. The solution gives access 
to all local quantities of the original lattice Hamiltonian 
Tis-f + SW-s-i, both static and dynamic. However, in this 
more complicated situation we were not able to exploit 
the cavity method as in Ref. |T2] and relate local quan¬ 


tities to non-local ones so to directly calculate (cJ ct c jo .). 
We thence opted for an indirect route and calculated the 
free-energy of the impurity model F lmp (h) as function of 
the source field h, shown in Fig. [S] for different tempera¬ 
tures below (metal) and above (insulator) the Mott tran¬ 
sition at U = 2.7 D. However, even though the change 
that occurs at the finite temperature Mott transition is 
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observable, the dependence on h is not regular enough 
to state unquestionably whether or not the impurity sus¬ 
ceptibility diverges in the limit h —» 0. Therefore we 
cannot give any definite answer to the question about 
the finiteness of the gauge-field expectation value. 



FIG. 8: (Color online) Impurity free energy for the model 
with the source term of Eq. ( |41| ). We plot the results vs. h 
for a range of temperatures, both in the metallic and in the 
insulating phase. 


IV. CONCLUSION 

We have shown that the Mott metal-to-insulator tran¬ 
sition in infinitely coordinated lattices can be faithfully 
reproduced by a model of fermions coupled in a Zi gauge- 
invariant manner to Ising spins in a transverse field. 
We numerically exactly solved such a model by dynami¬ 
cal mean field theory using the numerical renormaliza¬ 


tion group as impurity solver, and found that gauge- 
symmetry breaking spontaneously occurs at zero tem¬ 
perature and corresponds to the metal phase in the orig¬ 
inal electronic model. The zero-temperature metal-to- 
insulator transition translates in this spin-fermion rep¬ 
resentation to the recovery of the Z- 2 . gauge symmetry, 
thus endowing the Mott transition of a genuine order pa¬ 
rameter. Since we also found that gauge invariance is 
recovered at any however small finite temperature, we 
suspect that the Ising magnetization order parameter is 
closely related to the Fermi surface discontinuity in the 
quasiparticle distribution, which also disappears at any 
non-zero temperature. While the exact phase diagram 
agrees qualitatively well with that obtained by mean- 
field theory, the dynamical properties differ substantially. 
Indeed, contrary to mean-field, we have found that the 
slave-spin spectral function is gapless everywhere in the 
metallic phase, where it shows a linear dependence at low 
frequency, and becomes gapped only in the Mott insula¬ 
tor. This is the case both at zero and at finite temper¬ 
ature, even though in the latter case the transition, for 
T < T c ~ 0.026 D , or the crossover, above T c , are not 
anymore related to the Ising spontaneous magnetization, 
which is zero for any T / 0. Our attempts to uncover 
the behavior of the gauge field, which is hidden in the 
theory, in order to better characterize the finite temper¬ 
ature phase transition/crossover were unfortunately not 
conclusive, even though the dynamical behavior of the 
auxiliary Ising fields notably changes across it. 
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